Cold matter trapping via slowly rotating helical potential. 
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We consider the cold bosonic ensemble trapped by a helical interference pattern in the optical loop 
scheme. This rotating helical potential is produced by the two slightly detuned counter propagating 
Laguerre-Gaussian laser beams with counter directed orbital angular momenta ±lh. The detuning 
5uj may occur due to rotational Doppler effect. The superfluid hydrodynamics is analysed for the 
large number of trapped atoms in Thomas-Fermi approximation. For the highly elongated trap the 
Gross-Pitaevskii equation is solved in a slowly varying envelope approximation. The speed of axial 
translation and angular momenta of interacting atomic cloud are evaluated. In the T — >■ limit 
the angular momentum of the helical cloud is expected to be zero while toroidal trapping geometry 
leads to 2lh angular momentum per trapped atom. 

PACS numbers: 37.10.Gh 42.50.Tx 67.85. Hj 42.65.Hw 



I. INTRODUCTION 

The hydrodynamics of the sufficiently cold (T ^ 
1Q~^K) bosonic ensemble trapped by optical potential 
V^(r, t) [l|-y| is described by the Gross-Pitaevskii equa- 
tion (GPE) [3] for macroscopic wavefunction '^{f,t) : 



9* h^ 
at 2m 



Airh' 



1^1'*, (1) 



where m is the mass of atom, as is the two body s-wave 
scattering length. The negative a^ reduces the energy 
of ensemble and causes the mutual attraction of atoms. 
This results in formation of bright solitons in ID and 
collapse in higher dimensions. On the contrary in repul- 
sive BECs atoms repel each other and dark solitons or 
vortices are formed. In the periodic potential gratings 

i: 

F(r,i) -/(z,r,i) ~ exp [~r^] cos[ 5 ujt-{kf+kb)z], (2) 

where r is a distance from propagation axis z, 6uj = 
c ■ (kf — fcf,) — iOf — ujb is a frequency difference of the 
counter-propagating z-paraxial laser beams with the op- 
posite wave vectors |fc(/,b)| ~ A:(/^{,), the many-body non- 
linearity leads to the nonlinear tunneling, self-trapping 
and other quantum interference phenomena [6[. For ac- 
celerated ID optical gratings when 5uj = Slo ■ t, i.e. 
in the non inertial reference frames [51] the interacting 
bosons demonstrate Bloch oscillations and Landau-Zener 
tunneling. In rotating ID lattices the repulsive ensem- 
ble may have stable ground states and vortex soliton 
states m. For the attractive ensemble in 2D rotating 
lattices '§] the localized stable solitons exist in the cer- 
tain range of angular velocities. Inside the parabolic well 
with trapping frequency u!± and rotation frequency fl the 
nonlinear localized matter waves appear when n ~ uj± 
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at lowest Landau level [y] which appears due to the 
fictionary magnetic field induced by trap rotation 110!]. 
The goal of the present work is to study bosonic en- 
semble in rotating reference frame by a virtue of rotat- 
ing optical dipole trap whose rotation is due to the fre- 
quency detuning Slo. The proposed trap configuration 
is composed of the two counter propagating optical vor- 
tices. The shape of the interference pattern is defined by 
the mutual orientation of their orbital angular momenta 
±m (OAM). When OAMs are co-directed the A/2 spaced 
toroidal traps are formed by Laguerre-Gaussian(LG) [ll[ 
or Bessel vortices [12]. This geometry had been used for 
the single atom trapping and detection [13] and for the 
analysis of persistent condensate flows in LG beam waist 
14]. When OAMs are counter directed due to the phase- 
conjugation of the backward reflected LG beam the ID 
sinusoidal intensity grating is transformed into the t ruly 
3D helicoidal (figHlsecII) grating I{z, r, 9, t) [llt[ll[l5|. 
This grating experimentally observed for (5a; = [l6|: 

I{z, r, 0, t)=r'^^^^ exp [-r^] cos[ Sut - {kf + kb)z + 2£0], 

(3) 
where z,r,9 are cylindrical coordinates, must rotate with 
angular frequency H. = 9 = 5uj/2i when Soj ^ 0. 

Noteworthy the similar interferometric configuration 
with counter directed spin optical angular momenta (cir- 
cular polarizations) which is used for the sub-Doppler po- 
larization gradient cooling 17|. The interference pattern 
forms the static potential gratings due to Zeeman shift 
between ground-state sublevels. This causes the Sizyphus 
deceleration of the atomic beam in optical molasses. 

We consider two possible regimes of the cold ensemble 
trapping when kinetic energy is small compared to the 
interaction and trapping energies. One Thomas-Fermi 
solution is obtained as a balance of the red-detuned op- 
tical attractive potential and self-defocusing due to pos- 
itive as- This solution has atomic density p^ perfectly 
collocated with rotating optical helix I{z, r, 9, t). In this 
case the rotating potential imposes rotation to superfiuid. 
The other TFA solution has nonrotating density " funnel" 
profile Pf^^ relevant to thermodynamic limit when num- 
ber of atoms iV — ^ o3. 



Next the analytical solution of GPE for nonzero kinetic 
energy '^hei is described. It is valid when optical trap- 
ping and interaction energies are comparable and sub- 
tract each other in GPE. The linear momentum < P^ > 
and angular momentum < L^ > of this helical atomic 
cloud are evaluated. The Landau criterion \V\ > e{p)/\p\ 
for the appearance of elementary excitations and super- 
fluidity breakup [l^ is discussed for helical geometry. 



II. TWISTED WAVETRAINS 

It is well known that interference of a two counter 
propagating waves with a different frequencies a;/,aj(, 
roduces a running sinusoidal roll intensity grating 
[l9| . For the equal wave amplitudes |E/|,|Eb| and the 
phase difference (j) the distribution of the light intensity 
I{z, r, 9, t) has the following form: 

I{z, r, 9, t) - 2|E(/,b) 1^ [1 + cos[ Stjjt - (fc/ + fcb)z + ]] 

provided that a visibility of pattern is good enough 
|E(/)| ^ |E(b)| [2fl], where zr is Rayleigh range, Dq 
is a beam waist radius, the self-similar variable x = 
{iOf — LOb)t — {kf -\- kb)z + (j)] is responsible for the trans- 
lation of the interference pattern along z axis with the 
group velocity T4 — (uf — a;b)/(fc/ + fcb). The transversal 
(in the plane {r,0)) confinement of the light amplitudes 
E^,Eb is typical for the zeroth-order Gaussian beams. 
The roll interference pattern evolves into the sequence 
of the equidistantly spaced rotationally invariant in 9 
ellipsoids centered at the propagation axis z llj. For 



the higher-order propagation modes namely Laguerre- 
Gaussian beams (LG) with azimuthal quantum number 
i and orbital angular momentum ih per photon |2l| : 



E(^^b)(z,r,0,i) 



E(/^b) exp[i(-a;(/_b)t ± fc(/^b)z) ± ie9] 



[l+iz/zR] 



{r/D.r 



exp 



Da\l+izlzR) 



,(5) 



the interference pattern is different for LG reflected 
from conventional mirror and phase conjugating mirror. 
Backward reflection from conventional spherical mirror 
changes the topological charge of the LG [ll|, exactly 
in the same way as circular polarization of light changes 
from left to right and vice- versa in reflection [l^. The 
intensity I^^^ (z, r, ^, t) vanishes on the beam axis thus in- 
terference pattern transforms into the sequence of the 
equidistant rotationally invariant toroids separated by 
A/2 interval: 



I.A^,r,9,t) 
exp' 



An[l + cos[5u ■t-{kf + fcb)z]](r/Do)'l^l 
2r2 , . 2|E(y,b)|22(l^l+i) 



Dn 



'■{i + zyzir 



-, ^n 



TTe\a 



(6) 




FIG. 1: The isosurface of the optical intensity /^^ and the 
Thomas-Fermi density ph of the cold atomic cloud in a he- 
lical optical dipole trap (|7ll9p ( the scales are in ^m, but 
longitudinal modulation of A/2 is enlarged). The spatial mod- 
ulation is induced by the interference of counter-propagating 
LG beams with the opposite angular momenta. The pattern 
rotates with angular frequency 5(jj/2£ as a "solid-body" . Mag- 
netic field B adjusts the scattering length Ua to balance the 
attractive optical potential by many body defocusing. 



where £„ is the dielectric permittivity of vacuum. The 
reflection from phase-conjugating mirror (PCM) does not 
change the topological charge of LG and the interference 
pattern is twisted [ll|, [ig] : 



/,„(z,r,0,i) = [l + 

An ■ {r/Dof' 



cos[ Suj ■ t - {kf + kb)z + 2£ 9]] x 
2r2 



exp 



D^{l + z^/zl) 



(7) 



The intensity also vanishes at LG axis z as r^'^', while a 
self-similar argument: 



X = [{ujf- ujb)t - {kf + kb)z + 219], 



(8) 







keeps the maxima of intensity at the 2i collocated he- 
lices separated from each other by A/2 interval (fig. [Ij. 
The azimuthal term 2£ 9 appears due to phase conju- 
gation Eb '-^ E/*. Thus we have the following strict 
correspondence between the roll interference pattern ([2]) 
and the helical interference pattern Q: the frequency 
difference Slu = ujf — ujb is the cause of the translation 
of rolls with group velocity T4 — {tof — ujb)/{kf + kb) of 
the wavetrain, produced by a sum of the two counter- 
propagating beams (E/ -I- Eb) @. The 5(jJ is responsi- 
ble also for the rotation of helices with angular velocity 
9 = (5a;/2€. The rotation is the cause of the pitch of 
helical interference maxima along z-axis. Consequently 
there exists a perfect mechanical analogy between the 
solid body rotation of the helix described by formula ([7]) 
and an Archimedean screw. Namely the positive Suj cor- 
responds to the counter-clockwise rotation and this pro- 
vides the pitch in positive z direction for right helices. 
On the other hand the negative 5u) means clockwise rota- 
tion. In this case {5uj < Q) the positive translation speed 



in z direction takes place for the Ze/i-handed helices. Ev- 
idently the change of the topological charge i changes the 
direction of helix translation V^ due to alternation of the 
helix hand to the opposite one for a given 8io. 

This mechanical analogy is useful for the analysis of 
the cold atoms motion in the helical trap. The velocity 
vector of the condensate fragment trapped and perfectly 
collocated with intensity maxima has two components 
y = Fz + V6((fig. [1]). The axial component is a group 
velocity |T4| = (w/ — ujb)/{kf + k^) of the wavetrain, 
while the azimuthal component [Vel = {ujf — ujf,) ■ Dq is 
of kinematic nature. 

The electromagnetic orbital angular momentum inside 
LG beam waist volume V = ttDqZ]^ located near z = 
plane within Rayleigh range \z\ < zn is the expectation 
value of the angular momentum operator L^ — —ih [ fx 



<^->(/.b) = <*(/,6)l^-l*(/,6) > 



= 2eoJ {E+,^^y{-th[rx V]^(+_,))dV= 2eo> 
J{E+^^^ri-th-^E+,^) rdr-dddz c, ±lh 



hf,b}V 



(9) 



Hf,b) = eoC|E(/,f,)P is the light intensity, *^^,,j 



where /( 

y/'2eg Et,^^Az,r,9,t) are the macroscopic wavefunctions 
of a single photon inside a forward or backward beam 
([5|) with the winding number (., E~^ means the positive 
frequency components in E{r,t) spectrum. The square 



modulus 1^ 



if, by 



is a probability density of the photon 



detection which is proportional to the energy density of 
classical wave [2J]. In this particular paraxial case the 
spin-orbit coupling [25| is small enough and the angu- 
lar momentum of the photon is exactly decoupled to the 
spin and the orbital component: J = S + L. The lin- 
ear momentum expectation values for the forward and 
backward LG are as follows: 



<Pz> 



=<* 



(^, 



if,b)' 



if.by 



if,b) 



IPI* 



if,b) 



>= 2eo X 



±hk 



(f,b) 



'if,b) 



V 



hjjj 



U.b)^ 



(10) 



The ratio of the angular and linear momenta is Lz/ Pz ^ 
ic/ujrfi,\ [25]. The angular and linear momenta of the 
composite wavetrains ([6]) and ([7]) (fig. [T]) are: 



<Lz> 



^ ^if+b) 



^< ¥f + ■^t'lU'^'f + -^t' > 



{ih±eh) 



kf,b)V 



<Pz> 



(f+by- 



c^{hkf-hh)^^i^, (11) 

where the upper sign in ± corresponds to reflection from 
conventional mirror while the bottom sign stands for the 



reflection from PC-mirror with alternation of the photon 
angular momentum. 

The simplest conceivable configuration of linear PCM 
composed of plane mirrors, wavefront curvature compen- 
sating lenses and beam splitters is described in [2j, [23] . 
The goal of proposed setup is to counter-direct the split- 
ted LG beams. The four reflections are the necessary 
minimum. The case is that the incidence angles much 
above 45 degrees will distort both polarization and spa- 
tial structure of LG. In addition a small sliding of a beam 
along the reflecting surface 26] may occur when mir- 
ror is tilted with respect to the LG propagation axis. 
The evident physical restriction on this loop setup is 
to keep the path difference AL of the counter directed 
LG smaller than coherence length of trapping laser field 
(AL << c-Tcoh) [23]. The small frequency shift 5u] « 
27r • 10~^^^^'rad/ sec required to cause the helix rotation 
[ll | might be induced by a frequency ramp [5] or via ro- 
tational Doppler shift which appears due to rotation of 
the half- wavelength plate [2^ or Dove prism [23] . 



III. COLD ENSEMBLE DENSITY AND 
VELOCITY FIELD IN THE HELICAL TRAP 

Consider a bosonic cloud prepared in an elongated trap 
[a] and suddenly released afterwards. The well elaborated 
experimental procedure is to impose a periodic optical 
potential to study the Bloch oscillations, macroscopic 
Landau-Zener tunneling and Josephson effects [j, [a] . In 
our case the imposed optical potential is a helical one: 

VoA^,r,e,t)^-^^P^I^^{z,r,e,t), 



2e„c 



a{(jj) = 67reo c^ 



r/c^o' 



{ujI-u^-z{u^Iu^^)T) 



, (12) 



where a(w) is the polarizability of atom, which is real, 
i.e. a{(jj) « iieja] at large detunings from resonance 
Lo — ujf), V = e^ujQ^ /6TreQmi,c^ is classical damping rate 
via radiative energy loss, We is electron mass [3[. The 
GPE for wavefunction Vl' of ensemble confined by Vopt is: 



ih 



(9*(f,i) 



--—A* + Vopt 
2m 



m 

(13) 

with as{B) = a&g(l + Ab/{\B\ — Bp)) magnetic field 
dependent s — wave scattering length, where a^g is back- 
ground value of Os, Bp and A^ are the Feshbach mag- 
netic induction and resonance width respectively. 

Consider the sufficiently large number of trapped 
atoms {N = 10^~^^). Then quantum pressure 
term following from the uncertainty principle is small 
{h'^A'^ /2m ^ 0) compared to the optical trapping and 
interaction terms [4i] : 



Ei, 



N ■ abg 



Eki 



Oho 



(14) 



where Oho — yhjrrujjz is the " harmonic oscillator width" 
Q. The ratio Eint/Ekm (fT4|) is rauch more than unity for 
the most of the near infrared lasers A = 0.8 — 1.5 fim^ un- 
der the standard focusing requirement Dq ^ 10 — lOOfim 
and when as{B) is tuned in the range = 1 — lOOnm via 
Feshbach resonance. For the mid-infrared trapping at 
CO2 lasing wavelength A = 10.6/iTO the kinetic energy 
term Ekin = 2h^/mX^ is about 100 times smaller. The 
optical dipole trapping energy Edipoie is: 



Edipole — N ■ Vopt — 



2eoc 



=Aw2 2 



p- E 

Auj 



~aE^N 



eoc 



^f,b 



2 

-Wo 



(15) 



where Itw is the optical intensity of trapping beams [3| . 
Hence one might expect that a Thomas-Fermi approxi- 
mation (TFA) is adequate in our case when interaction 
is repulsive a^ > Q. 

When inhomogeneous ensemble is considered in the 
thermodynamic limit {N — >■ cx)) the local density ap- 
proximation for chemical potential is used fj,{r) = 
Miocii [p(^)] + ^optir) [1| . Let us consider the following 
TFA wavefunction in the vicinity of the LG-beam waist 
(i.e. within Rayleigh range \z\ < zn): 

^^„„ = $(r)exp [~^-^-^ + i^rf sm{Sujt-2k^z + 2ee)], 

(16) 
where the local r-dependent chemical potential fi(r) is: 

/x(r) = 47rn2a,$(r)Vm; lua = 2eoC|E(^,f,)|2; 

(17) 



^rf = exp(-2rVi?o^) • (r/i^o)^'^' • "^'"^'''" 



2e,-|C ■ h ■ Su! 

The wavefunction (ITB)) is normalizable and fits the 
GPE by substitution. This TFA density of atoms 
p^^Jz,r,e,t) = |*/„„|2 ^ expi-2{r/Dof)/ir/D„)^\'\ is 
{z,t)— independent "funnel" collocated with the optical 
helix I^^{z,r,9,t). The phase modulation of ^/,„i has a 
maximum near the density maximum p,^{z,r,9,t)^ si- 
nusoidal dependence on azimuthal angle 9 and decreases 
down to zero on LG axis and outside the LG waist. Note- 
worthy ^f„„ is a multiply valued function of 9. Hence 
this solution is of restricted interest. It may be used for 
evaluations of the thermodynamical_ parameters of the 
cold ensemble with |^/imp density [4|. 

The other solution for dilute Bose gas (e.g. N ^ 10^) 
is obtained in the TFA for r-independent chemical po- 
tential /.t(r) = const. The mean field wavefunction here 
^h is a sum of the two phase-conjugated vortices with 
the opposite angular momenta ±M: 

*, (z, r, 9, t) = *,(z, r, 9, t) + *_,(z, r, 9, t) ^ 



^±dz^O} 



1 + zV4 



exp 



exp[ ^ — h ikfz + U9] 

(1 + iz/zB) 



exp[ 



Do'{1+iz/zr) 
iii.t 



ikbz — i£9] 



(1 -t- Iz/zr) 



,(18) 



where the difference of the partial chemical poten- 
tials {p^ — A*j,)) associated with the each of "counter- 
propagating" wavefunctions "i'^^'i'^e is adjusted to the 
frequency difference of counter-propagating optical fields 
{fij: — p^)/h — 5 uj = uif — uJh- The substitution of this 
TFA wavefunction into GPE gives the following link for 
parameters: 

a(a.)/,, (0)-[l -f- co^jSut -t- 2kz ± 219)] ^ 2[,[ 

^^ o 71~, — 9T^2\ l'/-^0j ^ 



exp 



2e„c- 5 (1 + ^2/4) 
2r2 

D^il + z^zl) 



Ph{z,r,9,t), (19) 



where k — kf = k/,, p = fi w /i^^ is a constant 
(z,r, 6',i-independent) value of chemical potential, g — 
4Trh^as{B)/m is the interaction parameter. This is the 
quasiclassical restriction imposed on a homogeneity of 
chemical potential of the system in an external field Vopt 
[J] . In accordance to this solution the density of the cold 
atomic ensemble ph{z,r,9,t) is perfectly correlated with 
the rotating optical helix potential I^^{z,r,9,t) as de- 
picted at figHJ The density ph rotates as a " solid body" . 
The speed of the axial translation is Vz — X ■ Suj/Att. 

Apart from TF approximation the truly exact solution 
with nonzero kinetic energy exists '^ex = 'i'f + 'i'b- It 
is also a superposition of the two counter propagating 
paraxial matter waves ^/ and ^b. Let us reduce GPE 
(|13p to the paraxial form relevant to highly elongated 
geometry, as a helical one in our case: 



ih- 



a* 



(f,b) 



dt 



2m " ' 2m oz 



+Vopt *(/,fc) + — fc2(/,b)*(/,6) + 5l*/ + *6| '*(/,fc)(20) 

The perfect mutual cancellation of trapping and interac- 
tion terms: 

- Voptiz, r, 9, t) *(;^b) = 5l*/ + *fc| '*(/,fc), (21) 

occurs when nonlinear defocusing due to positive scatter- 
ing length as is compensated by attraction to intensity 
maxima caused by red detuning. Taking again the ex- 
act helical solution of (|20l) as a superposition of the two 
counter-propagating vortices: 

^^ Jz, r, 0, i) = */(z, r, 0, t) -f *b(z, r, 0, i) = 

*/-exp[-^ + ikfz] + *b-exp[-^ - ikbz], (22) 



and under the natural assumptions: 

k,f,b) d^U:b)ldz » d^^u.Jdz^ 



(23) 



two following equations for counter propagating and 
counter rotating matter waves "^ f and ^j, are valid: 



(24) 



which have the vortex solutions with charge i for the 
^^ t,(z, r, ff) with initial condition at z = equals ^o- 



*o-(r/£'o)'^'-exp 



I- 



(/,6) 



Da^(\+izlzii)) 



±iie 



(1 + iz/zp;) 

(25) 
The issues of dynamical stability of (P^ with respect to 
small perturbations and thermodynamical stability from 
the point of view of least energy arguments deserve a 
further careful analysis and will be published elsewhere. 
For example transformation to the reference frame rotat- 
ing with angular velocity Jl = 5ujI2(. synchronously with 
trapping helix leads to equation [J, |^, [Ifll : 



(9* ft2 



V, 



opt 



exp( 



-2r2 



Dl{l + z-^/k-^D^) 






where stationary solutions for the diluted (/^ = const) 
ensemble ^ = $(z, r, 6) exp{—iiit/h) are given by: 

ft^ - o (9<i> 

A*$ = -—A'i> + Vopt{z, r, 9) $+g|$| H+Qih—. (27) 



IV. MACROSCOPIC OBSERVABLES 

The helical solutions composed of the counter prop- 
agating free space LG wavefunctions apparently fit the 
continuity equation and have realistic field of velocities. 
In addition to the ensemble density p(r, t) obtained above 
in TFA the structure of velocity field V{'r, t) is a conse- 
quence of the complex geometry of the helical wavetrain 
which requires the perfect adjustment of the phase fronts 
achieved by phase-conjugation of colliding vortices. The 
one possible application of the (|18[ W^ might be in using 
them as variational anzatz for emulation of GPE [30j. 
Nevertheless the explicit form of solutions '^h (|18I22[) of- 
fers a possibility to evaluate the macroscopic observables 
of the trapped ensemble. Helical wavetrain has nonzero 
momentum P, (fTT 



< P^ >h^< ■^h\-^h—\^h >= Nh{kf -kh)^h dv- 

(JZ Jy 

exp{-r^)r^^^^[kf - kt + kfexp{ix) - kbexp{-ix)], (28) 
and easily calculated angular momentum L^ ^: 

< L, >h^< ^h\ - *^^l*'» >= ^^^(1 T 1) ^ 

h I dT/exp(-r2)r2|^l[l=Fl + exp(ix)Texp(-ix)], (29) 
Jv 

due to the apparent identity J^ sin(x)(i0 = 0. The up- 
per =p sign in (j29p corresponds to counter directed angu- 
lar momenta and helical interference pattern ([7]), while 



bottom =F sign corresponds to the toroidal optical inter- 
ference pattern ^. 

The same expectation values < Pz >ex—Nh{kf — fcf,) 
and < Lz >ex= N£h{l =F 1) has exact wavefunction 
^c^ (|22p . Quantum mechanically this happens because 
the wavefunction in both cases is a superposition of the 
two partial matter waves ^/ and ^h in ((22|) (or ^^ and 
^~i in (fT8|) ) having opposite and quantized (i.e. equal to 
±ih) mutually subtracted angular momenta. This means 
also that helical optical wavetrain ([7]) contains 0AM of 
exactly x £h per photon, while toroidal wavetrain ([6]) 
contains 2£H per photon PTj) as shown in sec.IL 

This feature looks seemingly counter intuitively from 
the point of view of classical hydrodynamics, but the 
similar results on vanishing of the moment of inertia 
for purely superfluid ensemble were summarized in [4|. 
Namely the density of atomic ensemble ph — \'^h,ex\'^ ro- 
tates as a solid body and one might expect that ph to have 
classically the angular momentum Ldass = Izz ■ Su}/2£, 
where Izz is the moment of inertia of the helical wire lo- 
cated in LG beam waist with the density profile ph [Sli]: 

Izz =- f Nm\^h^ex\''r^dV ^ Nm f |*,,,e.| V^drd^dz c^ 
- Nm {1 + cos{Sujt + {kf + kb)z + 2£0)) ■ dOdz ■ 



exp(-rVi?o')r- 



2+2|f| 



rdr - NmDo^Zr. (30) 



Nevertheless due to the quantization of the angular 
momentum in free space, the oppositely directed angu- 
lar momenta cancel each other completely, because they 
have integer opposite values of ±£S [S^] ■ In T — > limit 
[4| the net angular momentum of the helical wavetrains 
p8l[22)) is zero because the superfiuid component remains 
only. On the contrary, the linear momentum is not quan- 
tized in free space and this leads to nonzero net linear 
momentum Pz of the ensembles p8l22p , regardless to the 
mutual orientation of their OAM's. The net momentum 
Pz is small because of the smallness of the group velocity 
of the helical wavetrain Vz — duj/(kf + kt). For exam- 
ple when frequency splitting 6uj is induced by rotational 
Doppler effect [23, [28| , the speed of the axial translation 
of the helical density profiles (|18I22I) Vz ^ is several p,m 
per second (several rotations of helix per second) . 

The else interesting physical consequences relevant to 
experiments with trapped quantum gas traps may be for- 
mulated from the point of view of the Landau criterion 
\V\ > e{p}/\p\ for the appearance of elementary excita- 
tions (rotons) and superfluidity destruction, where e{p) 
is the energy - momentum dispersion relation. Following 
to 18, 22] consider the flow of quantum gas in a narrow 
helical channel with velocity V — Vz + Ve- In the rest 
frame the momentum of excitation p must be opposite to 
the velocity of superfluid V, because of the least energy 
constraint imposed upon excitation e{p)+p-V < 0. Thus 
e(p) — |p|-|y| < and |y| > e{p)/\p\. Because in our case 
the only significant component of VisVg = SujDo/2£ the 



excitations with monientum p appear when: 

5uJ,r^tDo>2i■e{p)/\p\. 



(31) 



The experimentahy controllable detuning 5uj of counter 
propagating waves w/ and Wf, by rotational Doppler ef- 
fect [23] which leads to the change the angular velocity 
of helix rotation makes possible to determine the critical 
velocity of superfluid, defined by contact point of roton 
minimum of e(p) with the line |p|Ve. The turbulent ex- 
citations (rotons) are assumed to appear due to ejection 
of superfiuid across the trapping potential barrier owing 
to centrifugal force, rather than because of the roughness 
or the channel end llM wA ■ 



V. CONCLUSION 

The flow of the degenerate quantum gas in helical 
trap had been studied analytically in the framework of 
the Gross-Pitaevskii equation. The necessary conditions 
were formulated for the appearance of the helical Bose- 
Einstein condensate flows due to the Thomas-Fermi bal- 
ance of the self-defocusing of condensate with positive 
scattering length a^ and "red" detuned optical dipole 
potential. The minimal achievable ensemble tempera- 
ture might be approximately evaluated as a recoil one 



Trecoii = 4 • h^/{2mX'^ ■ ks) [i3|- The possible exper- 
imental implementation of helical trapping is a sudden 
switching on of the helical potential after the conden- 
sate release from elongated optical trap in a way similar 
to switching of accelerated grating in Ref. [5] . 

The peculiarities of the cooling mechanisms in this heli- 
cal configuration were not considered in the current work. 
But the helical interference pattern (see figlj) geome- 
try might reveal the new features of the well elaborated 
mechanisms as a the Doppler cooling [33], polarization 
gradient cooling ^17] or velocity selective population trap- 
ping [SJl- The newly found loo p a nd helical features of 
the optical speckle patterns [3 a . |36| are also a promising 
trapping opportunities which may enlighten the features 
of the Anderson localization of cold atoms in ID and 3D 
speckle patterns |37| . 

Noteworthy the similar helical geometry of the collid- 
ing LG optical vortices of picosecond duration with oppo- 
site angular momenta proposed recently for the plasma 
currents excitation via ponderomotive force [23|. As the 
plasma vortices are the sources of the axial magnetic 
fields, the superfluid motion in helical trapping environ- 
ment (jl3|) is to be associated with a so-called artificial 
magnetic fields [3, [l3l ■ 
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